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We report theoretical calculations of high harmonic generation (HHG) by intense infrared lasers 
in atomic and molecular targets taking into account the macroscopic propagation of both funda- 
mental and harmonic fields. On the examples of Ar and N2, we demonstrate that these ab tnitto 
calculations are capable of accurately reproducing available experimental results with isotropic and 
aligned target media. We further present detailed analysis of HHG intensity and phase, under var- 
ious experimental conditions, in particular, as the wavelength of the driving laser changes. Most 
importantly, our results strongly support the factorization of HHG at the macroscopic level into a 
product of a returning electron wave packet and the photorecombination transition dipole, under 
typical experimental conditions. This implies that the single- atom/molecule structure information 
can be retrieved from experimentally measured HHG spectra. 
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I. INTRODUCTION 

High harmonic generation (HHG) is an extreme nonlin- 
ear optical process in vifhich an intense ultrafast infrared 
light is efficiently converted to an ultrafast coherent ex- 
treme ultraviolet (XUV) or soft X-ray light. In the last 
two decades, HHG has been widely studied for its poten- 
tial as a short-wavelength light source, either in the form 
of a useful bright, coherent tabletop light down to the 
water- window region (280-540 eV) [l,l2|, or the produc- 
tion of ultrashort light pulses such as single attosecond 
pulses and attosecond pulse trains [1, 01 ■ Recently, HHG 
itself has also been shown to have the potential to im- 
age molecular structure with sub- Angstrom precision in 
space and sub- femtosecond precision in time [5-8]. The 
basic principle of harmonic emission in a gas medium is 
well understood qualitatively. When an atom or molecule 
is exposed to an intense laser field, first, at a certain 
time an electron wave packet tunnels through the barrier 
formed by the combined atomic and laser fields. Next, 
it propagates in the laser field and may be driven back 
to recollide with the target ion. High harmonics are gen- 
erated when the returning electrons recombine with the 
ion, and convert the e nerg y gained in the laser field to 
high-energy photons [9l.ll0| Since the laser field interacts 
with a macroscopic medium, and the harmonics from all 
atoms or molecules are generated coherently, a full de- 
scription of the experimentally observed HHG spectra 
requires the treatment of the nonlinear propagation of 
the fundamental laser beam together with the harmon- 
ics in the medium. Thus the study of HHG consists of 
two parts. The first one is the calculation of the induced 
dipole by each atom or molecule in the laser field. This 
is to be carried out quantum mechanically by solving the 
time-dependent Schrodinger equation (TDSE) or equiv- 
alents accurately. The second one is to solve the non- 
linear propagation of the fundamental laser pulse and 
the harmonic fields in the medium by using Maxwell's 
equations. As in all nonlinear processes, an efficient har- 



monic generation requires good phase matching from all 
the elementary induced dipoles, which in turn depends 
on laser's properties like intensity, pulse duration, pulse 
shape, and target properties such as gas pressure, posi- 
tion of gas jet with respect to the laser focus, in addition 
to the linear and nonlinear responses of gas in the light 
fields. Clearly, a full understanding of HHG cannot be 
reached until all of these effects are properly described 
theoretically. This is especially important if one is to use 
HHG to image the structure of molecules. Since HHG 
spectra are sensitive to the detailed experimental condi- 
tions which usually cannot be accurately determined in 
a given experiment, how to extract useful quantitative 
structure information of individual target molecules in 
the gas medium is clearly a challenge. 

The most accurate way to obtain the induced dipole 
of an atom or molecule in a laser field is to solve the 
TDSE numerically. Since this approach is quite time 
consuming and the calculations have to be carried out 
for hundreds of laser peak intensities in order to describe 
the nonuniform laser distributions inside a focused laser 
beam, this is rarely done in existing studies including 
macroscopic propagation effect of HHG Instead, 
the much simpler strong field approximation (SFA), or 
the so-called Lewenstein model [l^l, is often used to cal- 
culate the single atom response. Despite of this limi- 
tation, the temporal and spatial properties of HHG ob- 
served experimentally have been reasonably understood 
from such SEA-based calculations. On the other hand, 
in a few examples, macroscopic HHG spectra obtained 
using TSDE-calculated induced dipoles do show signif- 
icant quantitative discrepancy compared to SEA-based 
calculations [H, [13] , and such studies have been limited 
to a few atomic gases only. To image the structure of 
molecules from the experimental HHG spectra, one first 
needs to be able to describe HHG spectra from molecules 
including macroscopic effects. 

In this article, we demonstrate an accurate and efficient 
method for calculating the HHG spectra from a macro- 
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scopic atomic or molecular gaseous medium. The method 
is based on our recently developed quantitative rescat- 
tering (QRS) theory |15l - [l7l | which allows us to calculate 
the induced dipole of an atom or molecule in a laser field 
with accuracy comparable to those obtained from solving 
TDSE, yet with computing time comparable to those by 
using the SFA. More importantly, according to QRS, one 
can express the complex induced dipole moment as the 
product of a complex returning wave packet and a com- 
plex photorecombination (PR) transition dipole moment, 
where the former depends on the properties of the laser 
and the latter is solely the property of the target. In fact, 
the elementary PR transition dipole moment is identical 
to the laser-free elementary photoionization (PI) transi- 
tion dipole moment which has been well studied in the 
last few decades. Using the QRS, we further show that 
the complex returning wave packet can be obtained from 
the SFA. The validity of the QRS, at the single atom 
or single molecule level, has been carefully calibrated 
against TDSE results for one-electron model atoms, and 
against experimental HHG spectra from molecules jT5| . 
Clearly such comparison is incomplete without consid- 
ering the macroscopic propagation effects. In an earlier 
paper, based on the laser induced dipoles calculated us- 
ing the QRS, Jin et al. [HI studied the macroscopic 
propagation effects of the HHG of rare gases theoreti- 
cally for the situation where the laser intensity and the 
gas pressure are small. Under this limit, the fundamental 
laser field is assumed not modified during the propaga- 
tion. It was shown that the macroscopic HHG spectra 
after propagation can be expressed as the product of a 
"macroscopic wave packet" and the same single-atom PR 
transition dipole. This result is very important since it 
enables us to extract target structure from the exper- 
imentally measured HHG spectra, thus paving the way 
for using infrared lasers for time-resolved imaging of tran- 
sient molecules. In the present article, we extend the 
work of Jin et al. [l^ to higher laser intensities and 
gas pressures at which the nonlinear propagation of the 
fundamental field is considered. We then examine the 
theoretically simulated HHG spectra of Ar and compare 
them directly with experimental data. We further ex- 
tend the method to include molecular targets, which are 
aligned or isotropically distributed. 

In Section II, we summarize the method and the es- 
sential equations for describing the macroscopic propaga- 
tion, and the calculation of single-atom or single-molecule 
induced dipoles in the QRS model. We also stress that 
the HHG spectra should be calculated for each specific 
experimental condition. Within the QRS, we can define 
a "macroscopic wave packet" (MWP) which will refiect 
the effect of lasers and the consequence of propagation in 
the medium. In Section HI, the results are shown and an- 
alyzed. First, we consider HHG spectra of Ar generated 
by 1200 nm lasers, and show that the experimental HHG 
spectra from 30-90 eV can be accurately reproduced the- 
oretically based on the QRS, but not on the commonly 
used SFA. In fact, the HHG spectra depend on the posi- 



tion where and how the XUV light is measured. We show 
how the two-dimensional HHG spectra, their global be- 
havior and individual single harmonics in the far field, 
depend on the gas pressure and the pulse length. We 
also study the spatial distributions of individual harmon- 
ics in the near field and in the far field. The phase and 
the amplitude of the harmonics after propagation are also 
analyzed, for harmonics along the propagation axis and 
off the axis. The harmonics are found to be always var- 
ied with these parameters and the different experimental 
conditions. However, we find that all the differences can 
be attributed to the different MWP's. Thus the depen- 
dence of MWP on the gas-jet position with respect to the 
laser focus, the degree of phase matching with respect 
to the gas pressure for individual harmonics are inves- 
tigated. Since phase-matching condition is also depen- 
dent on the wavelength of the laser used, we investigate 
how macroscopic HHG scales with the laser wavelength. 
The QRS has been used to obtain induced dipoles from 
molecules, so we extend the present work to molecular 
targets. We consider the isotropic and partially aligned 
molecules, and demonstrate that the experimental HHG 
spectra of N2 molecules from recent measurements and 
the present calculations are in good agreement. In Sec- 
tion IV, we summarize and discuss future perspective be- 
fore concluding this paper. 



II. THEORETICAL METHOD 

A. Propagation of the fundamental field 

In an ionizing gas, the propagation of a fundamen- 
tal driving laser is affected by refraction, nonlinear self- 
focusing, ionization, and plasma defocusing. The pulse 
evolution in such a medium is described by a three- 
dimensional (3-D) Maxwell's wave equation |19l - [2l| : 



V2^i(r,z,i) 



1 d^Ei{r,z,t) 



^(1 



dt 



(1) 



where Ei{r, z,t) is the transverse electric field of the fun- 
damental laser pulse with frequency wq. In cylindrical 
coordinates, = V^^ -I- d^/dz^, where z is the axial 
propagation direction. The effective refractive index of 
the gas medium can be written as 



77cff(r, z, t) = rioir, z, t) + ri2l{r, z, t) 



ujl{r,z,t) 



2^0 



(2) 



The first term 770 — 1 + (5i — takes into account of 
refraction {5i) and absorption (/3i) effects of the neu- 
tral atoms, the second term accounts for the optical Kerr 
nonlinearity which depends on laser intensity I{t), and 
the third term is from free electrons which contains the 
plasma frequency Wp = [e^no(t)/(£oWc)]^^^, where me 
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and e are the mass and charge of an electron, respec- 
tively, and ne{t) is the density of free electrons. The ab- 
sorption term J^hsi t) due to the ionization of the medium 
is expressed as [ij, [23 | 



Jahs{t) 



\EiitW ' 



(3) 



where jit) is the ionization rate, and Ip is the ionization 
potential. This term is usually small under the conditions 
for harmonic generation [l^ . 

The absorption effect (/3i) on the fundamental laser 
field caused by neutral atoms is in general small, so it is 
neglected. We only keep the real terms in the refractive 
index i]cs, and Eq. ([1]) can be written as 



dt 



+ -^Eiir,z,t)-2^{5i+T]2l)Ei{r,z,t). 



(4) 



By going to a moving coordinate frame (z' — z and t' = 
t — z/ c) and employing the paraxial approximation (i.e., 
neglecting d'^Ei/ dz'"^), we obtain [23j 



ViSi(r,z',0 



2d^Ei{r, z',t') 
c dz'dt' 



dt' 



-'-^Ei{r,z',t') - 2:^^(6, +r,2l)E,{r,z',t'). (5) 
c 

The temporal derivative in Eq. ([5]) can be eliminated by 
a Fourier transform, yielding the equation 



VlE,ir,z',uj) 



2iuj dEi{r, z',uj) 
c dz' 



= G(r,z',c.), (6) 



where 



and 



Ei{r,z',oj)^F[E,{r,z',t% 



(7) 



G{r,z', cj) =F{fio 



9 Jabs (r, Z',t') 

dt' c2 



Eiir,z',t') 



-2'^[6^+mHr,z',t')mr,z',t')}, (8) 

where F is the Fourier transform operator acting on the 
temporal coordinate. 

The plasma frequency ujp{r, z' ,t') is determined by the 
free-electron density nc{t'), and nc{t') can be calculated 
as following 

no(r, = no{l - cxp[- / 7(r, z', T)dT]}, (9) 

J —oo 

where no is the neutral atom density, and 7(r, z',t) is 
the ionization rate calculated from Ammosov-Delone- 
Krainov (ADK) theory [23 - [26j . The refraction coefficient 
Si , depending on the pressure and temperature of the gas 



medium, is obtained from the Sellmeier equation [27|, |28j . 
The second-order refractive index 772, also depending on 
pressure of the gas medium, can be calculated through 
third-order susceptibility x^^\ which can be measured 
from experiments ^29, 30] . Note that the relationship be- 
tween r]2 and x^^'' in Koga et al. [3l| differs from that 
in Boyd [ssj since the latter is derived by using time- 
averaged intensity of the optical field. The fundamental 
laser field is assmned to be Gaussian both in space and 
in time at the entrance of a gas jet (z' = Zin), and the 
pressure is assumed constant within the gas jet. 



B. Propagation of the harmonic field 

The 3-D propagation equation of the harmonic field is 
described by [3,[2l[3| 



V^E^{r,z,t) 



1 d'^Ei,{r,z,t) d'^Pir,z,t) 



(10) 



where P{r,z,t) is the polarization depending upon the 
applied optical field Ei{r, z,t). In this equation, the free- 
electron dispersion is neglected because the frequencies 
of high harmonics are much higher than the plasma fre- 
quency. Again going to a moving coordinate frame and 
using the paraxial approximation, Eq. (jlO|) becomes 



Vi£;h(r,z',t') 



2 d^Ei,{r,z',t') 
c dz'dt' 



d^ Pjr, z',t') 



(11) 



We eliminate the temporal derivative by a Fourier trans- 
form, obtaining the equation 



Vi^h(r, z',uj) 



2iuj dE\^{r, z' , lu) 



-uj^HQp{r, z',uj), 



where 



and 



^h(r, z',cj)^F[Ei,{r, z',t')] 



P{r, z',Lj) = F[P{r, z',t')]. 



(12) 



(13) 



(14) 



The source term on the right-hand side of Eq. de- 
scribes the response of the medium to the laser field and 
includes both linear and nonlinear terms. It is conve- 
nient to separate the polarization field into linear and 
nonlinear components P{r,z',uj) = x'"^K^)Ehii~, z' ,uj) + 
Pn\{r, z' ,uj), where the linear susceptibility x'^^H^) in- 
cludes both linear dispersion and absorption through its 
real and imaginary parts, respectively. The nonlinear 
polarization term Pni{r, z', lo) can be expressed as 

Pniir, z\io) ^ F{[nn - n,{r, z\t')]D{r, z\t')}, (15) 

where nc{r,z',t') is calculated from Eq. and 
D{r,z',t') is the single- atom induced dipole moment 
caused by the fundamental driving laser field. 
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The refractive index n{uj) = \/l + x^^K^)/^o [32 
related to atomic scattering factors by 



n{uj) 



(5h(w) - il3h{uj) 

-^noroA2(/i+i/2), (16) 

ZTT 



where rg is the classical electron radius, A is the wave- 
length, no is again the neutral atom density, and /i and 
/2 are atomic scattering factors which can be obtained 
from Refs. (s^, [sHi- Note that Sh{uj) and /3h(w) account 
for the dispersion and absorption of the medium on the 
harmonics, respectively. Finally Eq. (jl2p can be written 
as 



2iuj dEh{r, z',uj) 



(17) 



where the nonlinear polarization as the source of the har- 
monics is explicitly given. After the propagation in the 
medium, we obtain the near-field harmonics at the exit 
face of the gas jet {z' = Zput)- 

As presented in Ref. 18!] , when both the pressure and 
laser intensity are low, the fundamental field is not mod- 
ified through the medium. In other words, the source 
term in Eq. ([T|) can be taken as zero. Then the fun- 
damental laser field, assuming to be a Gaussian beam 
spatially, is given approximately in an analytical form. 
For the harmonic field, the dispersion and absorption ef- 
fects from the medium, which are explicitly expressed as 
a dispersion-absorption term in Eq. (|17l) are not included. 
These effects would become important if the gas pressure 
is high. For molecular targets, we will limit ourselves to 
experiments carried out under the conditions of low laser 
intensity and low gas pressure. Note that Eqs. ([6]) and 
([T7| are solved using a Crank-Nicholson routine for each 
value of uj. Typical parameters used in the calculations 
are 200~300 grid points along the radial direction and 
400 grid points along the longitudinal direction. 



C. Far-field harmonic emission 

Experimentally, harmonics are not measured at the 
exit face of a gas medium. They may go through a slit, an 
iris or a pinhole, or be reflected by a mirror before they 
reach the detector. The far-field harmonic emissions can 
be obtained from near-field harmonic emissions at the 
exit face of a gas medium through a Hankel transforma- 
tion (MSI 



El{r{,Z{,u;) 



. Eh{r,z',uj) km X 
I : — Jay- 1) 



Zi - 



X exp[ 



2{zf - z') 



Zf - 



rdr, 



(18) 



where Jq is the zero-order Bessel function, Z{ is the far- 
field position from the laser focus, r{ is the transverse 
coordinate in the far field, and the wave vector k is given 
hy k — UJ / c. 

Suppose the harmonics in the far field are collected 
from an extended area, the power spectrum of the macro- 
scopic harmonics is obtained by integrating harmonic 
yields over the area: 



5h(w) oc 



\El{x{,y{, Zf,uj)\'^dx{dyf. 



(19) 



where X{ and j/f are the Cartesian coordinates on the 
plane perpendicular to the propagation direction, and 



D. Quantitative rescattering (QRS) theory 

In this work, the single-atom (or single-molecule) in- 
duced dipole moment D{t') in Eq. ([T5l) is calculated 
quantum mechanically using the QRS theory. A detailed 
discussion of QRS for HHG from atoms or molecules is 
given in Ref. [l7] . We briefiy discuss the QRS theory for 
atoms and molecules separately in the following. 



1. Atomic target 

According to the QRS, the induced dipole moment 
D{uj) can be written as [SOj 



D{lu) = W{uj)d{Lu), 



(20) 



where d{uj) is the complex photorecombination (PR) 
transition dipole matrix element, and W{uj) is the com- 
plex microscopic wave packet. |W^(a;)p describes the flux 
of the returning electrons and is the property of the laser 
only. The QRS is a simple model that improves the SEA. 
It replaces the plane wave used in the SEA with accurate 
scattering wave in the calculation of PR transition dipole 
matrix elements, while the returning microscopic wave 
packet is the same as that in the SFA. Since the electron 
wave packet after tunneling but before the recombination 
is governed mostly by the laser field while the electron 
is far away from the target ion, and such interaction is 
fully described by the SFA, this explains why the elec- 
tron wave packet (its dependence on the momentum of 
returning electrons) derived from the SFA is accurate. 
In practical applications, the QRS obtains the induced 
dipole moment by 



(21) 



where both D^^-^{uj) and d^^^(w) are complex numbers, 
while d^^-^{uj) is either a pure real or pure imaginary 
number. Within the single active electron (SAE) approx- 
imation, we calculate ^'^^^(a;) using "exact" numerical 
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wave functions for the bound and continuum states. For 
Ar, we use the model potential given by Miiller [ioj . 



V{r) = -[1 + Ae-'- + (17 - ^)e-'^'']/ 



(22) 



with(A=5.4, C=3.682). In this model, spin-orbit interac- 
tion is neglected. The parameters have been chosen such 
that the minimum in the photoionization (or photore- 
combination) cross section is reproduced correctly. We 
comment that in principle the parameters in Eq. (|2ip can 
be generalized to many-electron wave functions if needed. 



2. Molecular target 

Within the QRS theory, the induced dipole moment 
D{uj, 9) for a fixed-in-space molecule is given explicitly 
by 



D(w, 9) ^ N{ef''^W{uj)d{uj, 9), 



(23) 



where N{9) is the alignment-dependent ionization prob- 
ability, W{uj) is the microscopic wave packet, and d{uj, 9) 
is the alignment-dependent transition dipole (complex in 
general) . Here 9 is angle between the molecular axis with 
respect to the laser's polarization. We limit ourselves 
here to linearly polarized lights, linear molecules and con- 
sider the parallel component of HHG with respect to the 
laser polarization only. Thus only the parallel component 
of the transition dipole d{uj,9) is needed in the calcula- 
tion. Note that W{uj) does not depend on the alignment 
angle 9, and it can be calculated formally as 



W{io) 



N{eY/^d{uj,9)' 



(24) 



Recall that W{ijj) can be obtained from SFA, where all 
the matrix elements above are calculated by replacing the 
continuum waves by plane waves. Since the wave packet 
W{lo) is independent of the alignment angle 0, it needs to 
be calculated only once for a given angle 9. In the QRS, 
the single-molecule induced dipole moment by the same 
laser is then obtained from Eq. (j24p by combining with 
the accurate d{ui,9) obtained from quantum chemistry 
code [4l| and with the tunneling ionization rate N{9) 
obtained from the MO- ADK theory [H, HI . Applica- 
tions of the QRS for HHG from single molecules have 
been investigated previously [13, IH, |4J] . 

Linear molecules can only be partially aligned when 
they are placed in a short laser field (pump laser) {is^. 
The intensity of the aligning laser is usually weak and not 
tightly focused such that it can be assumed to be con- 
stant within the gas medium. In other words, the degree 
of molecular alignment is not varied in the medium. In 
the following, the polarization of aligning laser is assumed 
to be parallel to the probe laser. The averaged induced 
dipole from the partially aligned molecules at each point 
in the gas medium is then obtained by coherently averag- 
ing induced dipole moments over the molecular angular 



distributions, i.e. 



D{uj,9)p{9) sin 9d9, 



(25) 



where p{9) is the angular (or alignment) distribution of 
the molecules. Similarly, the free electron density in 
Eq. is replaced by the averaged one: 



avg 



it') 



n^{t',9)p{9) ain9d9, 



(26) 



where nc{t' ,9) is the alignment-dependent ionization 
probability derived from Eq. For randomly dis- 

tributed molecules, p{9) is a constant. Once the aver- 
aged induced dipoles D^^s{uj) are obtained for a number 
of different laser intensities, they are then fed into the 
propagation equations for the harmonics. The propaga- 
tion is then carried out similar to that for atomic targets. 
We comment that in this model, dielectric properties of 
molecules due to non-isotropic distributions have been 
neglected. 



E. Macroscopic wave packet (MWP) 

As presented in Ref. [l^ [see their Eq. (25)], the 
macroscopic HHG spectrum in the near field for atomic 
targets can be expressed as 



(27) 



where W'{u!) is called a "macroscopic wave packet" 
(MWP), and d{uj) is the PR transition dipole moment. 
This relation still holds in the far field [see Eq. (IT5|) ]. 
since the PR transition dipole can be factorized out in 
Eq. The propagation of harmonics in free space to 

the far field would thus only modify the MWP. 

If molecules are only partially aligned, we calculate 
d(w) by coherently averaging the PR transition dipole 
weighted by the ionization probability of N{9) [25l. l42j: 



N{9y^^d{uj,9)p{9)siTi0d9. (28) 



From Ec^. (|27p . the target structure is reflected in the 
PR transition dipole, the propagation effect of the har- 
monics, in the meanwhile, is incorporated in the MWP. 
The two properties are well separated. The MWP rep- 
resents the cumulative effect of the returning electron 
wave packet (or microscopic wave packet) after propaga- 
tion in the medium and in the free space. The validity 
of Eq. ([TT)) forms the basis of extracting target molecular 
structure information from the experimentally measured 
HHG spectra. 
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FIG. 1: (Color online) Comparison of theoretical (QRS: dot- 
ted line; SFA: dot-dashed line) and experimental (solid line) 
HHG yields of Ar under a 1200 nm laser. Experimental data 
are from Ref. [4^. Laser parameters are given in the text. 



III. RESULTS AND DISCUSSION 

A. Macroscopic HHG spectra of Ar: theory vs 
experiment 

In Fig. [TJ we show the macroscopic HHG spectra gen- 
erated by a 1200 nm laser. Experimentally [4S], a 0.5 
mm-long gas jet was placed 3 mm after the laser focus 
(z=3 mm). A vertical slit with a diameter of 100 ^m 
is placed 24 cm after the gas jet. The beam waist at 
the laser focus is 47.5 /xm, and the laser pulse duration 
is ~40 fs. Laser peak intensity in the center of gas jet 
was 1.6x10^^ W/cm^, and the pressure of gas jet was 
estimated to be 28 Torr. In the simulation, the laser 
peak intensity and the pressure are adjusted until the 
best agreement with the experiment is reached visually. 
Using peak intensity of 1.5x10^'' W/cm^, we find the 
best agreement with cutoff position in the HHG spectra, 
and at pressure of 84 Torr, we find the best agreement in 
the widths of the harmonics [46| . 

We first calculate the single-atom response either by 
QRS or SFA. In SFA, Ar is treated as an effective "hy- 
drogenlike" atom where the nuclear charge is chosen such 
that its Is binding energy is the same as the binding en- 
ergy of Ar. The induced dipole moment is calculated by 
using the Lewenstein model [l3|. In the QRS, the "ex- 
act" transition dipole is calculated by using the model 
potential given by Miiller . The single-atom response 
is then fed into Eq. (flT)) . and the harmonic signals are 
collected in the far field in terms of Eq. ([T^ . 

In Fig. [U we can see very good agreement between 
QRS and experiment over the photon-energy region of 
30-90 eV. The "famous" Cooper minimum in Ar 47], is 
clearly seen in both experimental and theoretical spec- 
tra. The Cooper minimum in Ar has been reported in 
early HHG spectra generated by 800 nm lasers [H, 14^ , 
and it is seen more prominently using long-wavelength 
lasers Meanwhile, the propagated spectra obtained 
from SFA do not show Cooper minimum, nor does it 
reproduce the general spectral shape. Note that Cooper 
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FIG. 2: (Color online) Spatial distributions of harmonic emis- 
sion vs photon energy (normalized using on-axis intensity at 
77 eV) in the far field for the laser with different pulse du- 
rations and the gas jet with different pressures: (a) 24 fs, 84 
Torr, (b) 40 fs, 84 Torr, and (c) 40 fs, 28 Torr. The other 
laser parameters are not varied and are given in the text. 

minimum shall occur in the single-atom HHG spectra but 
it does not always appear in the macroscopic HHG spec- 
tra. As illustrated in recent simulations (46j , the position 
of Cooper minimum can change or even disappear under 
different experimental conditions. Such conclusions are 
consistent with experimental data where Cooper mini- 
mum may disappear in the HHG spectra by changing 
the gas pressure [4g| , or by changing the gas-jet position 
with respect to the laser focus [51|. In the following, we 
will show that these changes are due to variations in the 
MWP, and the separability of Eq. (gT]) is still valid. 

In Fig. [TJ there are still small discrepancies between 
the experimental data and the simulation by QRS de- 
spite of various attempts using somewhat different laser 
parameters. The harmonic width (or harmonic chirp) in 
the simulation is narrower than that in the experimental 
data. The harmonic width is mainly determined by laser 
intensity, pulse duration, and gas pressure [52h54| |. In 
the experiments, parameters like pressure of the gas jet, 
laser intensity and its spatial distribution cannot be mea- 
sured precisely. Other factors, like the use of the slit and 
the position of the detector can also influence the HHG 
spectra. All of these uncertainties can contribute to the 
discrepancy between the simulation and the measured 
HHG spectra. On the other hand, as seen from Fig. [1] 
the overall agreement between the experiment and the 
simulation over the 30-90 eV region is quite satisfactory. 

B. Harmonic chirp 

Although Fig. [T] shows the general global agreement 
between simulation and experimental data, it is pertinent 
to examine typical individual harmonics more carefully. 
How does the spectral width (or the harmonic chirp) vary 
as the pulse duration, gas pressure and laser intensity are 
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FIG. 3: (Color online) The dependence of spectral distribu- 
tion of individual harmonics on experimental conditions. The 
spectra are integrated over the vertical dimension as shown 
in Fig. [21 Experimental conditions are: 40 fs, 84 Torr (red, 
solid lines); 24 fs, 84 Torr (green, dashed lines); and 40 fs, 28 
Torr (blue, dot-dashed lines). Note that the peak position of 
the harmonic in each frame has been shifted to coincide for 
easy comparison. 



changed? In Fig. [T] the harmonics were taken for the gas 
jet placed after the laser focus, thus short trajectories 
were selected. In this case, the harmonic chirp is less 
dependent of the laser intensity, especially for a long- 
wavelength laser [s^ . We actually vary the laser intensity 
by 20%, the harmonic width almost does not change (not 
shown). We will concentrate on the effects of the pulse 
duration and gas pressure here only. 

In Figs. [2Ja)-[5l^c), we show the spatial distribution of 
harmonic emission in the far field (24 cm after the gas 
jet) by varying laser duration and gas pressure. All the 
other parameters are kept the same as in Fig.[TJ For each 
harmonic, the distribution on the vertical plane is shown. 
Integration of harmonic intensity over the vertical dimen- 
sion in Fig.[2lb) gives the simulated HHG spectra by QRS 
in Fig.[TJ A longer pulse duration and/or a lower pressure 
tend to generate sharper (narrower width) lower-order 
harmonics. A careful inspection reveals that the peak 
position of the harmonic actually blue-shifted from one 
frame to another. The shift is due to the change of the 
fundamental pulse as it propagates through the nonlin- 
ear medium In addition, the higher harmonics are 
less sharp, reflecting that the quality of phase matching 
varies with harmonic orders. 

In Figs. [3Ua)-|3td), we show the spectral distributions 
of harmonics H27 (27th harmonic), H37, II57, and H77, 
respectively, after they have been integrated along the 
vertical dimension. Here we examine the change of har- 
monic width as the pulse length and/or gas pressure are 
varied. Recall that the harmonic peaks are blue-shifted 
differently for different conditions. For easy comparison, 
the peak position is taken to be from the 40 fs, 84 Torr 
set (red, solid lines). The spectra from the other two sets 
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in the far field. Inset: Phases of H35-H65 in the region close 
to the propagation axis. 



are shifted to have the same peak position. From the 
figures, for H27 and H37, clearly the harmonic width in- 
creases with decreasing pulse duration. For a given pulse 
duration, the width increases with gas pressure. These 
figures also show that phase matching is not good for the 
higher harmonics, especially for the short-duration pulses 
where high-energy photons are emitted only from a few 
half-cycles. For these higher harmonics, narrower width 
seems to be obtained by raising the gas pressure. 

The harmonic chirp is a direct consequence of tem- 
poral variation of laser intensity, which can be mea- 
sured by XFROG (cross-correlation frequency resolved 
optical gating) [HI, 113]. It is determined by Auj{t) = 
—d^{u!,t)/dt. The harmonic phase ^{uj,t) is propor- 
tional to laser intensity I{t), with larger proportional 
constant for electrons taking long trajectories than for 
short trajectories [H^ [see their Eq. (1)]. For the focus- 
ing conditions in the present case, only short trajectories 
are selected and longer pulse leads to narrower harmonic 
width. In our model, Kerr effect on the fundamental 
field and plasma effect due to free electron are included 
in Eq. The ^{uj,t) is dependent on the gas pres- 

sure through ri2 and nc{t). It can be understood that the 
only variation of gas pressure could lead to the change 
of harmonic width. Similar analysis of the dependence 
of phase- matching on gas pressure can be found in Refs. 

MM- 



C. Harmonic divergence 

We next examine the harmonic emission in the near 
field and in the far field. The laser parameters are the 
same as those in the simulation in Fig. [T] Figs. HJa) 



8 



andlUJb) show that harmonic emissions in the near field 
are quite messy spatially ^Sj]. Because the harmonics 
are generated from a nonuniform Gaussian beam, phase- 
matching condition in the medium varies spatially. The 
radial variation of the phase ^{Lu,r) introduces a cur- 
vature of the phase front, which makes the harmonic 
emission divergent [58]. After further propagating to 
the far field, the harmonic emissions become regular, 
see Fig. SJ^c), where harmonics display Gaussian distri- 
butions centered at the propagation axis. In Fig. IDJc), 
only the short trajectories are selected. The harmonic 
emission is mainly on axis because of the small diver- 
gence. Since long- wavelength laser is used [6l|, the di- 
vergence in the region from H35 to H65 does not change 
too much. Similar study of the divergence of harmonics 
has been done for short-wavelength lasers [6^ . 

We also study the phase of harmonics vs the radial 
distance in Fig. U^d). For each harmonic, we note that 
the calculated phase grows quadratically with the radial 
distance and scales almost linearly with the frequency 
of the harmonics. The phases near the propagation axis 
vary for different harmonics, see the inset of Fig. Hfd). 
Below we show that the behaviors of intensity and phase 
of harmonics in the far field display good (laser-like) spa- 
tial Gaussian character. 

Recall that an incident Gaussian beam focused at z=0 
propagating along the z axis [l^, HI, is given by 



b + 2iz 



b + 2iz 



) = |i?(r,z)|e^^('-'^),(29) 



where 



\Eir,z)\= ,^f.°,^i/2 ^^P(~ fc2^_^4;,2 )' (30) 



(62-f 4z2)i/2 



and 



(/3(r,z) = -tan ( — ) + 



b ^2 _|. 4^2 ■ 



(31) 



Here £'o is the peak field at the focus, k is the wave vec- 
tor, b — 2t:wo/X is the confocal parameter, wq is the 
beam waist at the focus, and A is the wavelength. The 
intensity of each harmonic in Fig. UJc) follows the square 
of Eq. ([50)1 . The phase increases quadratically with r and 
linearly with the harmonic order, as seen in Fig.|3Kd), can 
be seen to follow the second term of Eq. ([3T|) . Near r=0, 
the phase from the first term of Eq. (I5T|) also contributes. 
This term (multiplied by harmonic order) gives a phase 
between -tt and tt for each harmonic, as seen in the inset 
of Fig. m^d) , in which the phase at r=0 is only defined 
within 27r. Since for different harmonics the confocal pa- 
rameters are probably changed either due to beam waist 
or wavelength, and the focused position may also change, 
we can only claim that each harmonic beam after propa- 
gation is close to a Gaussian beam qualitatively, but not 
necessarily quantitatively. Similar study of the phase of 
the harmonics in the near field was presented in Ref . [63| . 



D. Harmonic phase vs photon energy in the far 
field 



In the QRS model for single-atom response, the phase 
of each harmonic is the sum of the phase from the mi- 
croscopic wave packet and the phase from the transition 
dipole [see Eq. ([^0]) ]. For the harmonics calculated from 
the SFA, it is the same sum except that the phase of the 
transition dipole is either or vr. Thus the difference in 
the harmonic phase calculated from the QRS and from 
the SFA is given by the phase of the transition dipole 
calculated within the QRS (modulus tt). Is this relation 
still correct after the propagation, as implied by Eq. (|27l) ? 
In Figs. [5Ka)-[Sl[c), we show that phase differences (the 
squares) calculated by QRS and by SFA for far-field har- 
monics, and compare them with the phase (/'(w) of the 
transition dipole from the QRS (the solid line), at r=0 
mm, 0.5 mm, and 1.5 mm, respectively. The laser param- 
eters are the same as those in Fig. [T] The two are in good 
agreement. This agreement implies that the phase of the 
macroscopic wave packet W'{x, y, uj) obtained from QRS 
and from SFA remains identical after propagation at any 
points {x, y) on the plane. (In next subsection we will 
show that the magnitude of the MWP is also the same.) 

For the generation of attosecond pulses, the phase dif- 
ference (or atto-chirp) between consecutive harmonics is 
crucial [64j. The above results indicate that atto-chirp 
calculated using QRS and SFA differs only by the differ- 
ence of the phase 0(a;) of the transition dipole between 
two neighboring harmonics. This difference divided by 
2ojQ (in units of eV), or the derivative of 4>{lo) with respect 
to u) (in units of eV) is shown (dashed line) in Fig.[51[a). It 
is clear that correction to the atto-chirp calculated from 
SFA is small, except in the region near the Cooper min- 
imum where the phase 4>{uj) changes rapidly. This result 
is very significant since it explains why the generation of 
attosecond pulses can be understood mostly based on the 
SFA theory, even though it does not predict the harmonic 
spectra accurately. The phase of each harmonic is mostly 
determined by the phase of the returning electron wave 
packet, which has been accurately accounted for by the 
SFA, with very small corrections from the recombination 
process. This simplification explains why it is possible 
to study the generation of attosecond pulses in the last 
decade without a quantitatively accurate theory of HHG. 

To appreciate the complexities of harmonics, in 
Figs. [5l[d)-[5l[f) we show the far-field harmonic emissions 
at different radial distances r from the propagation axis. 
The harmonic emissions for different photon energies are 
comparable close to the axis, to provide a broad energy 
region from the on-axis area to synthesize attosecond 
pulses i65] . In the experiment this is accomplished by em- 
ploying an iris [M]- We note that at r=0 mm and r=0.5 
mm, the harmonic spectra resemble each other and the 
broad Cooper minimum can be easily identified. This is 
not the case if the spectra are taken at r=1.5 mm where 
the signals are much weaker and the Cooper minimum 
is no longer visible. As mentioned in Sec. Ill A, the 
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FIG. 5: (Color online) Phase difference (open squares) between far-field harmonics by QRS and by SFA, in comparison with 
the phase </>(aj) (solid lines) of "exact" transition dipole, at (a) r=0 mm, (b) r=0.5 mm, and (c) r=1.5 mm. The difference 
of atto-chip (the phase difference divided by the energy difference between two neighboring harmonics) between the QRS and 
the SFA is shown in (a) (dashed line). The far-field harmonic emissions at (d) r=0 mm, (e) r=0.5 mm, and (f) r=1.5 mm vs 
photon energy. Note that the phases in (a)-(c) are only taken at harmonic peaks, which are shown in (d)-(f), respectively. 



"disappearance" of Cooper minimum is attributed to the 
change of MWP, not due to the recombination process. 
Experimentally the "disappearance" of Cooper minimum 
in the HHG spectra of Ar has been shown using 800 nm 
pulses by changing the experimental conditions [H, [5l| . 



E. The dependence of MWP on experimental 
conditions 

In Fig. |6l[a) we show how the MWP (complex ampli- 
tude) changes with the gas-jet position calculated by the 
QRS, for HHG signal collected after the slit. For simplic- 
ity, only the envelope of the MWP is plotted. We keep 
laser intensity at 1.6x10^^ W/cm^ and gas pressure at 
56 Torr. The slit is placed 24 cm after the gas jet, and 
all other parameters are the same as given in Fig. [T] We 
also show the MWP calculated by SFA under exactly the 
same conditions, and the results are shown (in circles) in 
Fig. Ufa). The comparison shows that the MWPs (mag- 
nitude) from QRS and from SFA are the same for the 
same experimental condition, even though the MWP can 
change greatly depending on the gas-jet position. When 
the gas jet is placed "after" the laser focus (z=3 mm) 
the MWP is very flat, since good phase-matching is fa- 
vored for this arrangement as the single-atom harmonic 
phase is partially cancelled by the Gouy phase from the 
focused laser [gl, . If the gas jet is placed before the 
laser focus, we note that the MWP changes rapidly, es- 



pecially near photon energy around 50 eV. Such strong 
energy dependence can wash out the Cooper minimum 
in the HHG spectra j46| . 

In Fig.injb) we compare the MWPs as in Fig.lHJa) but 
without the slit. In such comparison, the harmonics from 
electrons following long and short trajectories are both 
collected. By comparing Figs, ^a) and IH^b), we note 
that the MWPs for gas jet at z—-3 mm and z—Q mm 
change significantly, but for z— 3 mm, the MWPs in both 
cases remain rather flat. This shows that a narrow slit in 
the far field can select the short trajectories effectively. 

In Fig. injc) we investigate how the MWP depends on 
the gas pressure for the focusing condition of z=3 mm. 
The MWP has been normalized by the ratio of the pres- 
sure. The three curves would be on top of each other 
if a complete phase-matching condition had been full- 
filled. The curve for higher pressure is slightly lower in- 
dicates that full phase matching is not reached, especially 
for the lower harmonics. With the increase of pressure, 
the MWP is much smoother vs energy. In fact, increas- 
ing the gas pressure tends to smooth out the harmonics. 
These results also indicate that the harmonic energy in- 
creases quadratically with the gas pressure, which is in 
agreement with measurements reported in Ref. 67[. In 



Fig.EJd) we display the fuh MWP vs the gas pressure. It 
shows that the harmonics exhibit better shaped peaks as 
the phase-matching condition is favored at higher pres- 
sure 
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FIG. 6: (Color online) (a) The magnitude (not the intensity) of macroscopic wave packets (MWP) at different gas-jet positions: 
z=-3 mm (before the focus, solid line), z=0 mm (at the focus, dashed line), and z=3 mm (after the focus, dot-dashed line) by 
QRS. The corresponding ones by SFA are also shown (open circles). The HHG signals are collected at the exit of the slit, (b) 
Same MWP as in (a), but the total HHG signals are collected without the slit. The curves in (a) and (b) have been normalized 
separately, (c) Dependence of MWP on gas pressure. Each curve has been normalized by the pressure, (d) The detailed 
structure of MWP in (c) from 55 to 85 eV. The curves have been shifted for easy visualization. 



F. Wavelength scaling of harmonic efficiency 

One of the main interests in the study of HHG is to 
produce bright tabletop XUV or soft X-ray sources, or 
intense attosecond pulses. Since the single-atom har- 
monic cutoff energy is proportional to the square of the 
wavelength of the driving laser, HHG generated by mid- 
infrared (MIR) lasers has been of great interest exper- 
imentally. While MIR lasers can efficiently reach high- 
energy photons, the yield is less favorable. It is of interest 
to study how the HHG yield scales with the laser wave- 
length. Within the single-atom response level, there have 
been a few theoretical calculations [68-73]. However, to 
compare with experimental HHG spectra, macroscopic 
propagation effect has to be included. A few investiga- 
tions on the wavelength scaling of HHG experimentally 
[sol IgtI Irsl - fTSj have been reported. However, theoretical 
analysis is still rather scarce. 

To study wavelength scaling of the HHG yields, one 
has to fix all other parameters that may affect the effi- 
ciency of HHG. One also has to decide if it is the total 
HHG yield or only the HHG yield within a given photon 
energy region. In single-atom simulations, the laser pa- 
rameters can be easily fixed. But this is not the case in 
experiments. Theoretical simulations including macro- 
scopic propagation effect are few [t^. Since the result- 
ing HHG spectra depend on so many other parameters. 



as we have demonstrated in the earlier subsections, any 
wavelength scaling laws derived are likely to depend on 
experimental parameters used. In spite of this limitation, 
it is still of interest to take a look at the wavelength scal- 
ing by using the present QRS model. For this purpose, 
we will define a parameter that describes the efficiency 
of harmonic generation. This is the ratio between the 
output energy (total harmonic energy) with respect to 
the input energy (fundamental laser energy) for different 
laser wavelengths. According to Ref. |2(|, the input en- 
ergy Epuiso can be related to the laser duration Tp and 
peak intensity Iq at the focus by iSpuisc — -^Triygrp if 
the laser beam has a Gaussian distribution in time and 
space. The output energy can be obtained by integrating 
the harmonic intensity over a photon-energy region: 

E'out = J J J \Ehix,y,uj)fdxdyduj. (32) 

In Fig.[71[a) we show the single-atom HHG spectra cal- 
culated for three wavelengths. Only the envelope of each 
spectra is shown. In the calculation, the laser intensity 
and duration are kept as l.GxlO^'* W/cm^ and 40 fs, re- 
spectively. In Fig. (Tljb) the HHG spectra obtained after 
including macroscopic propagation are shown. In the cal- 
culation, the beam waist at the focus is kept as 47.5 /im, 
a 0.5 mm- long gas jet is placed at 3 mm after the laser 
focus, and gas pressure is kept at 56 Torr. The yield of 
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FIG. 7: (Color online) (a) Single-atom HHG spectra, and macroscopic HHG spectra without (b) and with (c) the slit for 800 
nm (solid lines), 1200 nm (dashed lines) and 1600 nm (dot-dashed lines) lasers. The laser parameters are given in the text, (d) 
The wavelength dependence of the total integrated HHG yields above 20 eV. The integrated HHG yields in (a), (b), and (c) 



follow A-^■''±"■^ A-« "" ", and A-i«-2±0-2, respectively. 



each harmonic is obtained by integrating over the whole 
plane perpendicular to the propagation axis. In Fig.[7l^c), 
the HHG yields are recorded after they have passed a slit 
(the slit with a diameter of 100 fim is placed at 24 cm 
after the gas jet). From Figs[7Ub) and[7l^c), we calculate 
the HHG efficiencies per atom vs the wavelength. 

In Fig. [T^a), the ratio of input energy is 1:1:1 for 800 
nm, 1200 nm, and 1600 nm lasers. If we integrate the 
HHG yields above 20 cV as the output energy. The re- 
sulting energy follows ^-3-5±o.5 sJjq-^^q [j^ Y[g [7][d). If 
we integrate the HHG yields between 20 eV and 50 eV, 
which would give a scaling rule of A^^. In Tate et al. [68| . 
the laser intensity and the number of optical cycles were 
fixed for 800 nm and 2000 nm lasers. According to our 
approach, the ratio of input energy is 1:2.5 for the 800 
nm and 2000 nm lasers. And their scaling rules at con- 
stant intensity of A^^^^^-* would be modified as A^*^^^^^ 
at a constant input energy. 

When propagation effect is considered it is generally 
known [63, uM ^^^^ phase-matching condition is more 
difficult to meet for longer wavelength, thus the HHG 
efficiency decreases with increasing wavelength. Here 
we consider the total HHG yields for the lasers used in 
Fig. El^b) in which the gas jet is placed at z=3 mm. 
Since the laser intensity is fixed at the center of the 
thin gas jet, we calculate that the intensities at the 
laser focus are 1.78x10" W/cm^, 2.01x10^ W/cm^, 
and 2.33x10" W/cm^, for 800 nm, 1200 nm, and 1600 
nm lasers, respectively, thus the input energies have the 
ratios of 1:1.13:1.31. We find that HHG yields integrated 



from 20 eV up scale like ^-8-5±o.5^ g]-^Q-^j-^ Fig.[7l^d). 
If we only integrate the harmonics between 20-50 eV, 
then the scaling rule is A~-^°-^. 

Experimentally, Colosimo et al. [13] reported that the 
HHG yields between 35-50 eV for 2000 nm lasers are 
about 1000 times smaller than that for 800 nm lasers, for 
experimental conditions that were kept "as fixed as pos- 
sible" . This would give a A~^ dependence in this narrow 
energy region which is not too far off from our scaling of 
In addition, Shiner et al. ^&?\ reported a scaling 
rule of A^^ '^^^-^ for the HHG of Xe with a fixed laser in- 
tensity. By assuming perfect phase-matching condition 
for all the laser wavelengths used, they derived scaling 
law of A~^ '^^^-^ that was to be compared to the scal- 
ing law derived from the single-atom response. We can- 
not compare their results with our simulation. They also 
used a Bcssel beam (instead of Gaussian beam) in the ex- 
periment and the gas jet was located at the laser focus. 
Since the HHG yields depend on so many experimental 
parameters, it is clear that any simple scaling laws de- 
rived should be taken with caution. In Fig. [T^d), we also 
show the scaling law for the case where the HHG yields 
are collected after the slit. We integrate the HHG sig- 
nals above 20 eV and obtain the ;\-iO-2±o.2 s^-aijj^g^ jjj 

general, good phase-matching condition becomes more 
difficult to meet with increasing laser wavelength. Even 
if the gas jet is placed after the laser focus, the short tra- 
jectories are not selected efficiently for longer wavelength 
lasers. A slit is usually used to select contributions from 
short trajectories in the far field. By blocking out con- 
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tributions [see Fig. from the long trajectories the 

harmonic efficiency becomes worse. 

Based on the above analysis, the HHG yields for long- 
wave driving lasers under the same experimental condi- 
tions appear quite unfavorable. On the other hand, for 
practical purpose, experimentally high harmonics are to 
be gen erated with optimized conditions. In Colosimo et 
al. [50| it was reported that the HHG yields between 35- 
50 eV generated by using 2000 nm lasers can be as high 
as 50% of that from 800 nm lasers if the experimental 
conditions were optimized independently. Furthermore, 
Chen et al. [zl] demonstrated that it was possible to use 
much higher pressure to generate HHG for long wave- 
length lasers, thus achieving usable photon yields even 
in the water-window region. Clearly additional theoreti- 
cal analysis of macroscopic propagation effects on HHG 
for long-wavelength driving lasers under different exper- 
imental conditions is desirable. 



G. Macroscopic HHG spectra of N2 in an 800 nm 
laser 



The QRS has been used to calculate HHG spectra from 
single molecules. In order to compare with experimen- 
tal data, it is often assumed that HHG measured in the 
experiment are taken under the perfect phase-matching 
conditions. While such a model has been shown to be 
successful in in terp reting a number of experimental ob- 
servations [13, l43l . [4^ , it is still crucial to understand 
the effect of propagation on the HHG in the medium. In 
Figs.IDJa) andlSUb) we show the macroscopic HHG spec- 
tra generated by N2 molecules that are isotropically dis- 
tributed or partially aligned along the polarization axis 
of an 800 nm laser. The HHG spectra have been re- 
ported recently [zl] using 800 nm and 1200 nm lasers, 
and the results from 1200 nm laser have been recently 
analyzed f46j. To achieve good agreement with experi- 
ment in the cutoff region, we need to use peak intensity of 
l.%y.W^ W/cm^ instead of 2.3x10^'' W/cm^ in the ex- 
periment. We use the other parameters as close as those 
given in the experiment: pulse duration is ^ 32 fs, beam 
waist at the focus is ~ 40 /zm, and the slit with a diam- 
eter of 100 /im is placed at 24 cm after the gas jet. A 
1 mm-wide gas jet is located 3 mm after the laser focus. 
Using the QRS, we first obtain induced dipoles for fixed- 
in-space molecules for different laser intensities. They are 
then averaged coherently according to the alignment dis- 
tribution of molecules by Eq. (^5)) . The resulting induced 
dipoles are then fed into Eq. (|17p . In the experiment, the 
degree of alignment was estimated to be (cos^ 6')=0.6- 
0.65, we use an alignment distribution of cos* 9 in the 
simulation. Note that only HOMO is included in the 
calculation. This is adequate since contributions from 
HOMO-1 is important only for molecules that are nearly 
perpendicular to the polarization axis [3, [sO] • 

Figs. |51[a) and [SJb) show the good overall agreement 
between experiment and theory for both randomly dis- 



tributed and partially aligned N2. The experimental 
spectra reveal a weak minimum at 39±2 eV [79|] for both 
aligned and unaligned molecules. The theory also pre- 
dicts a minimum near 45 eV, and the position of mini- 
mum is not shifted from random to the partially aligned 



N2. For aligned N2, we note that McFarland et al. [81 1 
reported a weak minimum around 39 eV in the HHG 
spectra using an 800 nm laser. But the minimum was 
observed at about 45 eV in Torres et al. [Ill using an 
800 nm laser, and at about 38 eV using a 1300 nm laser 
for unaligned N2 (see their Fig. 3). For 1200 nm lasers, 
minimums in HHG spectra from the experiment [79] and 
the simulation |46j have also been reported. The exact 
location of the minimum is not always identical since it 
can be somewhat altered due to the energy dependence 
of the MWP [see Fig. [Hla)], which changes with laser 
parameters and experimental conditions. Despite of the 
difference in the positions of the minima between the 
simulation and experiment, we consider that the overall 
agreement in the two spectra is quite satisfactory. 

In Figs.[Sl[c)-[51[e) we show the phase difference between 
neighboring harmonics for randomly distributed and par- 
tially aligned N2 molecules in the far field, for different 
radial positions from the propagation axis. As mentioned 
in Sec. HI D, phase differences, which reveal emission 
times of harmonics [s^ and can be measured by RABITT 
(reconstruction of attosecond beating by interference of 
two-photon transition) [84, 85] technique, are significant 
for the generation of attosecond pulses [s^. Near the 
axis, we note that the phase difference is linearly changed 
vs photon energy in the plateau ending up by a sharp rise 
near about 48 eV in Fig.[51[c) and 45 eV in Fig.[5Kd), and 
beyond the cutoff it is almost a constant. The nearly 
constant phase difference was observed by Mairesse et 
al. [87] for atomic targets to optimize the conditions for 
attosecond pulse generation. Since the phase changes 
rapidly near the minimum of the harmonic spectra, this 
explains that the position of the HHG minimum in the 
integrated spectra can be easily changed, depending on 
how the integrated spectra are measured. We comment 
that the phase difference vs photon energy in Fig. [8][e) 
is less regular. These harmonics, taken at position away 
from the axis, have large contributions from long tra- 
jectories and they are not suitable for attosecond pulse 
generation. Finally, phase difference for randomly dis- 
tributed and partially aligned N2 agree well. The small 
variation comes from the slightly different derivative of 
the phase of averaged transition dipoles for random and 
partially aligned N2 [see Fig. |9Kb)]. 

To understand the results shown in Fig. [51 we show 
the MWP for the laser used, and photoionization (PI) 
transition dipole moment (parallel component only) for 
N2 molecules fixed in space, as well as the averaged PI 
transition dipoles defined by Eq. ([^ for randomly dis- 
tributed and partially aligned N2 molecules in Fig. [9l 
According to Eq. ([?7)) . macroscopic HHG yields can be 
expressed as the product of a MWP and an averaged 
transition dipole. The MWP (magnitude) is identical for 
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FIG. 8: (Color online) Comparison of theoretical (dashed lines) and experimental (solid lines) HHG yields from (a) random 
and (b) partially aligned N2 molecules in an 800 nm laser. Experimental data are from Ref. [t^. Phase difference between 
consecutive harmonics of random (solid circles) and partially aligned (open triangles) N2 molecules in the far field at different 
radial distances: (c) r=0 mm, (d) r=0.5 mm, and (e) r=1.5 mm. 
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FIG. 9; (Color online) (a) Macroscopic wave packet (MWP) 
for an 800 nm laser (solid line), and averaged photoioniza- 
tion (PI) transition dipole moments (parallel component) for 
randomly distributed (dashed line), partially aligned (dot- 
ted line), and perfectly aligned (^=0°, dot-dashed line) N2 
molecules, (b) The photon-energy dependence of the phase 
for averaged PI transition dipoles in (a). 



and the laser polarization axis (see Refs. [13, Sa])- An 
average over the angular distribution of the molecules 
washes out the minimum, except for a relatively faster 
drop of the transition dipole near 35-45 eV. The effect of 
angular average also washes out the rapid phase change 
in the transition dipole moment. For molecules fixed at 
6 = 0°, Fig. [9ljb) shows a rapid phase change of near tt 
in a very narrow energy region near 50 eV, i.e., at the 
position of the minimum in Fig. |9l[a). However, after 
the angular average, one can only see somewhat faster 
phase change at small photon energies. Since both the 
MWP and the transition dipole exhibit minor energy de- 
pendence, the actual minimum position of the HHG is 
difficult to locate accurately. With much better aligned 
molecules, the position of the minimum can probably 
be better determined. According to Eqs. (|27l) and (|28l) . 
the averaged PI transition dipole can be obtained from 
the experimental HHG spectra, and it could be used to 
retrieve the alignment-dependent ionization probability, 
N{9). This may provide another method to check the 
calculated N{9) 0. 



IV. SUMMARY AND OUTLOOK 



random or partially aligned N2 under the same laser and 
experimental condition. The transition dipole for fixed- 
in-spacc molecules shows a minimum at photon energy 
that depends on the angle between the molecular axis 



In this paper we have described a complete theory 
for high harmonic generation (HHG) in a macroscopic 
atomic or molecular medium. Our approach is based on 
the simultaneous solution of the coupled Maxwell's equa- 
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tions describing macroscopic propagation of both driving 
laser pulse and its high harmonic fields together with the 
microscopic induced dipoles. For the latter we use the 
recently developed quantitative rescattering (QRS) the- 
ory for a single-atom/molecule response. This scheme 
provides a simple and efficient method for calculating 
HHG from a macroscopic medium, which is otherwise 
formidable, especially in the case of molecular targets. 

Our results show quantitative good agreements with 
recent experimental HHG measurements [4^, [t^ either 
for Ar or N2 targets. For different laser and experimen- 
tal conditions, we present the detailed analysis of HHG 
intensity and phase. Since the calculation which includes 
macroscopic propagation for an isotropic or an aligned 
molecular target is quiet scare |!46], we hope that this 
paper will further stimulate the interest in establishing 
quantitative theory for HHG, which can compare directly 
with real experiments for partially aligned media. We 
note that in this paper the effect of absorption and dis- 
persion of high harmonics are neglected for molecular N2 , 
which should be adequate at low gas pressure. In gen- 



eral, absorption and dispersion are anisotropic in case of 
aligned molecules, which are not available generally in 
the literature. One of the most important results of this 
paper is that at the macroscopic level under typical ex- 
perimental conditions HHG spectra can be factorized as 
a product of a macroscopic wave packet and a photore- 
combination transition dipole. The latter is a property 
of the target only. This factorization provides a solid 
foundation for extracting the molecular structures from 
HHG spectra. As demonstrated recently \Jl, HHG has 
been used for ultrafast probes of excited molecules, our 
work also provides the needed theoretical basis for that. 
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